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Nonabelian gauge field theories coupled to fermions in some representation of the gauge group are asymptotically 
free as long as the number of fermion flavors is less than a certain number. Within this allowed range of fermion 
flavors, the theory is expected to be conflning in some range at the lower end and it is expected to be conformal 
at the higher end. Identiflcation of the critical number of fermion flavors that separate the confining region from 
the conformal region is a non-perturbative task that has recently received considerable attention within the lattice 
field theory community [I]-[T1]. Several issues need to be resolved before such an endeavor can make some physics 
conclusions. These include (a) How does one deal with a conformal theory on a lattice; (b) How does one conclude 
that a certain model exhibits features of near conformal features; (c) How does one compute the location of the 
infra-red fixed point in a lattice model. Since finite-volume effects need to be understood carefully and since one has 
to be close to the chiral limit to understand the above issues, the numerical simulations are inherently large scale in 
nature. 

An attractive alternative has been proposed to study large- iV gauge field theories coupled to adjoint fermions on a 
single-site lattice. It has been argued that reduction to a single-site lattice should hold in such theories [TS] and tested 
numerically using a variety of methods to see if one can reproduce the continuum infinite- volume theory by working on 
a single-site lattice [.16 - 28 . With the exception of [2Z) .28^, all attempts have considered the Eguchi-Kawai reduction 
and numerically argued that the single-site theory is in the correct continuum phase. Asymptotic freedom is maintained 
in these theories if the number of Dirac flavors is less than ^ . Two-loop perturbative beta function would suggest the 
existence of an infra-red fixed point if the the number of fermion flavors is greater than y|. With this in perspective, 
the single-site model with one massless adjoint overlap-Dirac fermion was extensively studied in [5S]. Numerical results 
suggest that the coupling runs much faster than what is predicted by continuum two-loop perturbation theory at the 
lattice couplings that were considered. In order to better understand the connection between single-site lattice models 
and infinite-volume continuum theories, we decided to revisit the problem of perturbation theory on the single-site 
lattice in this paper. 

We will consider the weak-coupling limit and the only parameters we will consider are the number of fermion flavors 
which we will extend to take on all real values in the range [O, ^] and the fermion mass. The main aim of this paper 
is to use a four-dimensional density function to answer two questions: 

1. What is the range of fermion flavors for which the single-site massless theory can be expected to reproduce the 
infinite-volume continuum theory? 

2. Can we reproduce the infinite-volume continuum theory with massive fermions? 

We will provide an answer to both these questions using Wilson fermions and overlap fermions. We will not consider 
the case of twisted reduction in this paper. 
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I. INTRODUCTION 
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II. THE SINGLE-SITE MODEL 



The single-site partition function for a SU(A^) gauge theory coupled to / flavors of fermions in the adjoint repre- 
sentation is 

Z = f U'^U.e^'^f^^ (1) 
with Haar measure dU. The Wilson gauge action is 

d 

Sg^bNj2 Tr [U^U^Ulut - l] . (2) 

The fcrmion action is 

S"/ = lndeti?„,o (3) 

with the subscript w for Wilson fermions and o for overlap fermions. The Hermitian Wilson Dirac operator for massive 
adjoint fermions is given by 

where mm is the bare Wilson fermion mass. The adjoint gauge fields are given by 

Af=^-Tr[T^U^T''Ul], (5) 
where T"^ , a — 1, . . . , {N'^ — 1) are traceless Hermitian matrices that generate the su{N) Lie algebra and satisfy 

Tr TT'' = 26"'' ; [T", T**] = ^ i/^^'T^ (6) 

c 

The Hermitian massive overlap Dirac operator is defined by 

Ho{mo) = ^ [(1 + mo) 75 + (1 - "^o) e [Hn,{m^)]] , (7) 

where vrio € [0, 1] is the bare overlap fermion mass and < is the irrelevant Wilson mass parameter. 
The total action depends on d SU(A^) matrices and the gauge transformation is 

U^^gU^g^. (8) 

Note that the eigenvalues of are gauge invariant. We cannot fix a gauge such that one of the 11^ = 1 since we are 
on a single-site lattice. The action has an additional U'^(l) symmetry given by 

C/^^e'"-f/^ (9) 

with < < 27r. Restricting to with integers < < N keeps it in SU(iV); otherwise we have trivially 
extended the SU(A^) theory to a V{N) theory. The four Polyakov loop operators, given by 

P„=Tr[/^, (10) 

are gauge invariant but not invariant under If the U''(l) symmetry is not broken, then the eigenvalues of all 
are uniformly distributed on the unit circle and = (the reverse statement is not necessarily true because the 
eigenvalues in different directions might be correlated). In the following, we set the number of Euclidean space-time 
dimensions c? to 4. 
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III. LEADING-ORDER PERTURBATION THEORY AND THE DENSITY FUNCTION 



The symmetry given by ([9| is spontaneously broken in the weak-couphng hmit if wc do not have adjoint fermions 
even when there are a finite number of flavors of fundamental fermions [29j . We want to study if this symmetry is 
spontaneously broken in the weak-coupling limit in the presence of adjoint fermions. We set 



and expand around V|^^ = 1 to compute observables in perturbation theory. 

The expression for the partition function at leading order at weak coupling is known |20| and is given by 



(11) 



5" = 5° + /5?: 



^Inp 



2^4sm 



(12) 



The fermionic contribution is 



5°,„ = 2 5^1n7;%( 



(13) 



where we have removed a ^-independent term that arises from the zero modes for massless fermions and assumed that 
pi-3 ^ if i 7^ j. The non-zero modes are given by 



sin^ [e^^-e^^) 



1 + mi 



1 - mj; TOn, -I- ^ 



(14) 



Owing to the symmetry given by (|9]) S'^ is invariant under — + ol^ for any choice of a^. 

As N ^ CO, we assume that we can define a joint distribution, p(0), in the following sense: At any finite N , for a 
fixed choice of i 



9^, i = 1, . . . , TV and = 1, . . . , 4, let 



(15) 



where 5 denotes the 27r-periodized delta function normalized to d6S{9) = 1. We can then rewrite S'^ in (12) as 



oO 



N'4d''ed^(bpie)Sgj{e~^)p{(b); 



P 



^4sh 



2 



Sf{6) = 21n7.u,^o(m^,o) ; 

.-X 2 

p 



Iwirriw) ^ im^ + -] +p; P = X! 



(16) 



Since there is a restriction in the sum that appears in (12) and (13), we have to evaluate the principal value of the 
integral appearing in (16) by excluding a small region around 6 = cj). The integral, J^, indicates the Cauchy Principal 
Value. Finally, we can write 



(17) 
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where by J [dp] we mean the integral over aU possible choices for djj^, j = 1, . . . , N and = 1 



We now assume that, as iV — > oo, the integral in (17) will be dominated by a single distribution p{6), maximi zin g 
S'''[/9]. We will only allow distributions that are non- negative everywhere with the normalization condition in (151 



Furthermore, we assume that the dominating distribution p{6) is smooth and finite for all 9 (in contrast to p defined 
in (151 for angle configurations at finite N). Since the singular nature of S{9) in (16 1 is only logarithmic^, the 
integrals are then finite even if we drop the principal- value restriction 6 ^ (p. Clearly, S*" in ( 16 ) is invariant under 
p{6) — p{6 + a) for any choice of a, corresponding to the invariance under ([9]). 
Owing to the periodic and symmetric nature of Sgj{9), it follows that 



Jo „ ^ 



(18) 
(19) 



Therefore, Fourier expanding 



> (27r)4 ^ 



results m 



with c_fc = cl, Co = 1 



(20) 



(21) 



provided p{9) is such that we can interchange the order of principal-value integration and sums over Fourier modes 
when we insert (20) in ([l6|. (If this is not the case, e.g. if p{9) is of the form (15), we expect the infinite sum in (21 ) 



to be diverging.) 

If all the eigenvalues, 



Xk = X 



ia) 



(22) 



for fc 7^ are smaller than zero, the constant mode, p{9) 



will dominate in the large- TV limit (i.e., — > for 



k ^ Q) and the single-site model will be in the correct continuum phase and possibly reproduce the infinite-volume 
continuum theory. In the next section, we will obtain the region in the (f^m^) plane for Wilson fcrmions and in the 
{f,mo,mw) space for overlap fermions where this is the case. Focusing on certain points in the allowed space we will 
compare the infinite- A'^ action from ( 16 ) with a numerically obtained maximum of the finite- A'^ action in (|12p to get 



a feel for the size of the finite- A'' effects. 



If some of the eigenvalues are larger than zero, then the action 5° in (161 will not be maximized by p{9) 



(2vr) 



and 



some Ck {k ^ 0) will be non-zero. Since the action in (21 1 is quadratic, the maximum will be obtained at the boundary 
of the domain of allowed values for the c^'s, which is determined by the condition p{9) > for all 9. Therefore, S[p] 
will be maximized by a p{9) which is zero at least at one point in the four-dimensional Brillouin zone. Due to the 
shift-invariance, there will then be a class of densities, related by p{9) — > p{9 + a) with arbitrary a, having identical 
maximum action resulting in a spontaneous breaking of the U''(l) symmetry in (joj). 



IV. INVESTIGATION OF THE ALLOWED REGIONS 



A. Overlap fermions 



We will start with the action for S**^ as given in ( 16 ) and find the eigenvalues Afc defined in ( 18 ) for all fc^ < 7 



(Xk is invariant under sign changes and permutations of the k^). In the following, we consider only fc 7^ 0. In 



order to compute the eigenvalues, we need to perform the integral in ( 19 1 numerically and we will do this using a 
four-dimensional uniform Riemann sum. 



^ A special case are massless fermions at / = 1/2, for which Sg{9) + fSf{9) is finite at = 0. 
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A sample plot is shown in Fig. [l] where we have computed the eigenvalues = A^^^ + fX^j/^ for massless overlap 
fermions with / = 1 and m-w = — 1- The results are obtained with equally spaced points in the four-dimensional 
integration space and we used M = 41 and M = 71 to show that we have reached the limit of the continuum integral. 
Since two eigenvalues are positive, (/ = 1, mo — 0, — —1) is not a point in the allowed region for overlap fermions. 
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FIG. 1: Eigenvalues Afc = A*^''' + /A^°' as a function of for the massless overlap Dirac operator with f — 1 and m„ = — 1 
obtained using numerical integration with equally spaced points in the four-dimensional integration space. 



As a second example, we set f — 2, keeping mo = and m^ — —1- In this case, we find all eigenvalues A^ to be 
negative, making this a point inside the allowed region. In Fig. [2j we have plotted ln(— Afc) as a function of to show 
that even in the log-scale we have a good estimate for the continuum integral. 



Numerically, we find that A^^^ > for all fc, which means that a point (/, TOq, m^) will be inside the allowed region 
(defined by A^ 



< for all k) iff 
< for aU k, 
(ii) />max,.{-Ai^VAi°^}. 



(i) \^^\mo,m. 



We observe that the eigenvalues A^^'' and A^"'' go to zero as fc -^• oo. As they approach zero from opposite sides 
(A'^^^ > 0, A^"-* < 0), we potentially have to consider all k in order to be able to determine the boundary of the allowed 
region in the {f,mo,my^) space. 



Let us first consider the case mo — 0. As fc — >■ oo, the integrals determining the eigenvalues in (19 1 are dominated 

g 



by (jji, € (0, 2T:/k^] since both Sg and So diverge as — >■ 0. Furthermore, expanding S{(j)) around (f> ~ 0, wc obtain 

^g(0) = -ln(02) + ... ; So{^) ^ 2ln{(f>^) + . . . , 



(23) 



indicating that — A). /A^ — >■ ^ as fc — >■ oo for all r7i„, < 0. Computing the eigenvalues numerically, we indeed find 

that — A^.^V'^i"^ rapidly converges to ^ for large fc (cf. Fig. [s] for an example). Therefore, for TOq = 0, the allowed 
region in the (m^,, /)-plane is determined by eigenvalues A^ with fc being small (cf. Fig. [s]). Considering only f < 



at fc = (2,2,2,2) for e [-1.21,- 
m^o e [-0.843,-0.780]. ^ [-1.21, 
(w„, /)-plane for mo = see Fig. |4] 



-1.15], at fc = (1,0,0,0) for m. 



we find numerically that the maximum — A^, /A^, , which leads to the boundary of the allowed region, is obtained 



e [-1.15,-0.843], and at fc = (2,0,0,0) for 
0.780] is not allowed. For a plot of the boundary of the allowed region in the 
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FIG. 2: Logarithm of the eigenvalues for the massless overlap Dirac operator with f — 2 and m^, = — 1 obtained using 
numerical integration with M** equally spaced points in the four-dimensional integration space. 
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FIG. 3: Plots of /o = — A^^'/A^°' (for < 9) at m-uj — —1 and different choices for nio. For rrio — 0, fo ^ 0.5 as — >■ oo; for 
all mo > 0, /o — >■ C)0 as fc^ — >■ oo. The boundary of the allowed region is determined by maxfc fo{k). 

For Too > 0, the divergence of So{(f>) at </> = is regulated and therefore x'"^\mo) / X^,^\mo = 0) — !• as fc — cx). 

Fig. [5] shows some examples for the dependence of a[.°''(too) on TOq and k at = — 1. Together with our results 

for the massless case, this immediately implies that — A^fV'^i"'' — 00 as fc — cx) (see Fig. [s] for numerical results). 
Therefore, it is necessary to keep TOq = in the weak-coupling limit. 
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FIG. 4: Boundary of the allowed region in the (m„,/)-plane for massless overlap fermions (mo — 0). The different lines 
show -A^^V^i°' for k = (2,2,2,2) (green, dashed), k = (1,0,0,0) (red, solid), and k = (2,0,0,0) (blue, dot-dashed). The 
intersection points are at (mu,,/) « (—0.843,2.30) and {rriu,,/) « (—1.15, 1.65). 

A'°V|A'°'(0)| 




FIG. 5: Plots of Xi°\mo)/\Xr{mo = 0)| at m„ = — 1 for modes with fc^j = fc V/i and different choices of k. The ratio goes to 
zero as fc — >■ oo for every rrio > 0, with the rate of convergence increasing with mo- 



B. Wilson fermions 

The scenario for Wilson fermions is very similar to the overlap case described above, with now playing the role 
of rrio and no additional irrelevant parameter. For m^; = 0, we have 15*^,(0) = 21n((/)^) + . . . for small 0, and therefore 

— A^^-'/A^^-' — > i as A; — > oo, while for > 0, we find —Al^V^fc"'' — J" oo as /c — ?> oo, which means that > is not 
allowed. Some numerical results are shown in Figs. [6] and [Tj which directly correspond to Figs. [3] and [s] for the overlap 
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case. For m^, = 0, the maximum — A^.^V'^i™'' is obtained at k = (1, 1, 1, 1), for which — A^^V-^i'"^ = 2.39 (cf. Fig. |6|. 
Therefore, only rriw — and / > 2.39 is allowed for Wilson fermions in the weak-coupling limit. 
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FIG. 6: Plots of fo = -A^^'/Ai"' (for < 9) for different choices of rriw For m^, = 0, /o — > 0.5 as — >• oo; for > 0, 
fo — >■ oo as fc^ — ^ oo. 
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FIG. 7: Plots of A^™' (rriu,)/! A^™' (m^, — 0)| for modes with fc^ = fc V/.1 and different choices of k. 



V. APPROACH TO THE INFINITE-iV LIMIT 



Given the leading-order partition function in (12 1, we view the action as a function of the AN angles subject to 
the condition that they belong to SU(A^) and perform a maximization of the action using the Hybrid Monte Carlo 
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algorithm as described in [50] . Instead of looking at the distribution of the angles which can look uniform due to the 
action being invariant under — >■ + for arbitrary integers fc^, we look at the action 5*" and compare to what 
one would get if we replace p{9) by the constant distribution T^prr in (16) which we will be N'^Xq. If the distributions 



at finite N given by ( |15[ ) for the maximum action configurations approach the uniform distribution as N ^ 00, we 
expect the action density, s° = to approach Aq. The plots in the top left panel of Fig. [s] are for two points in the 
allowed region and there is evidence for s° approaching Aq as TV — >■ 00. Contrary to this result, we see that s° does 
not approach Aq for two points outside the allowed region shown in the top right panel where we have only changed 
/ and kept mm = — 1 compared to the points shown in the top left panel. The bottom panel shows two other cases 
outside the allowed region where the parameters coincide with previous numerical work ^25) . 





FIG. 8: Approach of the lattice action density to the large-A'' limit. The left panel shows that we approach a uniform 
distribution in the allowed region. The right panel shows that we do not approach a uniform distribution in the region that is 
not allowed. The bottom panel shows two other cases where we do not approach a uniform distribution in the region that is 
not allowed (note that Aq < 0, i.e., lim]v->oo s" > A''). 



In order to provide f urth er support for the argument in the previous paragraph, we will start with the distribution 
at finite N as given in (15) for the maximum action configuration {0^} and compute all Fourier coefficients 



Ck 



(24) 



for all k with —3 < fc^ < 3. The results for several values of N for one point in the allowed region (/ — 2, nio — 0, 
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FIG. 9: The Fourier coefficients, \ck\, of p{9) as a function of fc = 
overlap fermions and the Wilson mass parameter set to mn, = — 1. 



'k^ for several values of TV with / = 2 flavors of massless 



TOtu = —1) that coincides with a point in the top left panel of Fig.jsjare shown in Fig. [oj Results for one point outside 
the allowed region (/ = 1, rrio = 0, = —4) that coincides with a point in the bottom panel of Fig. [S] are shown in 
Fig. [TOj We expect all the Fourier coefficients shown in Fig. |9]to approach zero and there is some evidence for this. 
We can imagine constructing a sequence of distributions for N = (with n = 2, 3, . . .) that approaches a uniform 
distribution in the large- A'^ limit by locating (5-functions on all sites of a four-dimensional periodic hypercubic lattice 
with lattice spacing 2tt /n. The corresponding Fourier coefficients would by 1 if all fc^ are multiplies of n and zero 
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FIG. 10: The Fourier coefficients, \ck\, of p{6) as a function of fc = Vk^ for several values of A'^ with / = 1 flavors of massless 
overlap fermions and the Wilson mass parameter set to mn, — —4. 



otherwise. It is therefore not surprising that we obtain non-zero Fourier coefficients Ck with k being of order in 
Fig. |9]even though we expect all coefficients to vanish in the large- limit. 

On the other hand, we expect some of the Fourier coefficients shown in Fig. [TO] to approach a non-zero limit and 
there is some evidence for this particularly when we look at the combined plot for N = 107, 157, 239 shown in the 
bottom right panel. 

If some of the Fourier coefficients shown in Fig. [TO] indeed approach a non-zero limit, the partial continuum action 



12 



density defined as 

s°= E ^fe^fc^fc' (25) 

|fe^|<3 

should not approach Xq. There is clear evidence for it when we look at the plots in the top right and bottom panels 
of Fig. |8] Note that the approach to TV — >■ oo is quite flat consistent with the convergence seen in the combined plot 
for N — 107, 157, 239 shown in the bottom right panel. Furthermore, the limit of and s° do not seem to coincide 
in the bottom panel of Fig. |8]for the case of / = 0.5 and — —3 suggesting that there are modes with k not in 
— 3 < fc^ < 3 that approach a non-zero limit at infinite N. 

Only i n th e limit of inflnite N, we are allowed to ignore the restriction of the principal value and sum the infinite 



series in (21) to obtain a finite action density provided the distribution has a smooth limit. The partial sum s", on 
the other hand, is finite at any N but will only agree with at infinite TV if all the coefficients not included in the 
sum approach zero excluding some accidental cancellation due to eigenvalues with different signs. If the distribution 
is uniform in the infinite- iV limit as is expected for points in the allowed region, we expect and to coincide with 
Aq. There is evidence for this in the top left panel of Fig. [Sj 



The computation of in ( 19 1 should exclude a small region of order e around </> = in order to properly account 



for the principal value required at finite N d ue t o the form of the distribution in ( 15 ). One has to tune e as a function 
of N and include a sum over all modes in J2l| (which will be finite) to match with at finite N. The difference 
between s° and at finite TV seen in Fig. plis a combination of two effects: not excluding a small region of order e 
and not including a sum over all modes. 

VI. DISCUSSION OF PREVIOUS NUMERICAL WORK 

Previous numerical work described in [20] only looked at one Fourier mode, namely, = (1, 0, 0, 0) and its permu- 
tations. We know from the analysis of the Fourier modes in this paper, that some coefficients could be accidentally 
small and it is necessary to look at several Fourier modes. The wilson mass parameter was set to rriw — —5 in the 
numerical analysis performed in [20 at finite lattice coupling and we know from the analysis performed here that this 
is not in the allowed region for any value of / in the weak-coupling limit. 

Numerical work in [5S] falls under a slightly different category. The running of the coupling studied in that paper at 
the range of lattice couplings did not agree with two loop perturbation theory. Since the theory studied used massless 
overlap fermions with / = 1 and rriw = —4, we know from the analysis performed here that we cannot obtain an 
infinite- volume continuum limit by going to — > oo (weak-coupling limit). The speculative part in [25j suggests the 
possibility of a continuum limit away from b ~ oo. If this is the case, then the analysis performed in this paper does 
not shed light into such a scenario. 

Numerical studies with two fiavors of Wilson fermions (both massless and massive) were carried out in p4] . From 
Fig. [6] in this paper, we know that we cannot obtain the infinite- volume continuum limit with one or two flavors of 
Wilson fermions. Evidence for being in the correct continuum phase was obtained by a study of operators of the form 
I TrC/^l, \Tr Uf^U,y \ and \Tr U^Ul\. These can be considered as special cases of 

Ck (b) = Tr C/^ U^^ U^' Ul' (26) 

which will tend to in ([24]) in the weak-coupling limit, 6 — > oo. The ordering of the operators will matter at finite 
lattice coupling if more than two of the fc^ are non-zero. We do not have the gauge field data for Wilson fermions 
used in [24^ but we do have them for massless overlap fermions at & = 0.65, / = 1 and m^, = —4 in j25j. The results 



for Cfe(6) are plotted in Fig. 11 and should be compared with Fig. 10 We see that b — 0.65 is far away from the 
weak-coupling limit consistent with the results in |25j . Furthermore, the coefficients with small k could be accidentally 
small. This suggests that the conclusions in [53] possibly result from being far away from the weak-coupling limit and 
not looking at a sufficient number of Fourier modes. 

In addition, we can also conclude that we cannot use a single-site model with heavy Wilson fermions and obtain 
the continuum limit of a pure gauge theory in contrast to the claims made in [26i . The eigenvalues of Wilson fermions 
are doubly degenerate j^D] but they do not come in pairs of opposite chirality. Therefore, it is not apriori clear how 
to deal with half a flavor of Wilson fermions making it essentially impossible to study the continuum limit of any 
infinite-volume theory using adjoint Wilson fermions on a single-site lattice. 

Since we cannot keep the bare mass finite and non-zero as we take the weak-coupling limit, the proposal in |21) to 
use single-site models with massive adjoint fermions in order to extract physics of pure gauge theories is ruled out. 
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FIG. 11: The Fourier coefficents, Ck{b), as a function of fc = Vk^ at finite lattice coupling, b = 0.65, with one flavor of massless 
overlap Dirac fermions and — —4. 

VII. FUTURE WORK 

The allowed regions plotted in Fig [4] provide for an interesting scenario when it comes to the usefulness of single- 
site theories to describe correct infinite-volume continuum physics. We cannot study theories with / < | unless we 
entertain the possibility that the continuum limit occurs away from b — oo but this would be a radical deviation 
from conventional wisdom. We can study / > 2 theories without a fermion mass. It is possible we can study / > 2 
theories with a fermion mass provided we can stay in the correct phase by taking the bare mass to zero as we go to 
the weak-coupling limit such that the physical mass is kept constant. Since this theory is expected to be conformal 
for massless fermions, it is not clear how a single-site theory will exhibit conformal behavior. This is certainly a case 
worth further investigation. 
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